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ABSTRACT 

Self-interactions of dark matter particles can potentially lead to an observable sep¬ 
aration between the dark matter halo and the stars of a galaxy moving through a 
region of large dark matter density. Such a separation has recently been observed in 
a galaxy falling into the core of the galaxy cluster Abell 3827. We estimated the DM 
self-interaction cross section needed to reproduce the observed effects and find that the 
sensitivity of Abell 3827 has been significantly overestimated in a previous study. Our 
corrected estimate is ct/todm ~ 3cm^ g“^ when self-interactions result in an effective 
drag force and cr/mDM ~ 1-5 cm^ g“^ for the case of contact interactions, in some 
tension with previous upper bounds. 

Key words: astroparticle physics - galaxies: clusters: individual: Abell 3827 - galax¬ 
ies: kinematics and dynamics - dark matter 


1 INTRODUCTION 


Recently Massey et al. (20151 have studied elliptical galax¬ 


ies falling into the core of the galaxy cluster Abell 3827, 
following previous such studies ( [Williams fc Saha|2011[ |Mo-| 
hammed et al. 20131. Using several strongly-lensed images 


of background objects, it is possible to reconstruct the po¬ 
sitions of the dark matter (DM) subhaloes of each of the 
galaxies. One of these is observed to be significantly sepa¬ 
rated from the galaxy’s stars by A = 1.62lQ'4g kpc. While 
noting that an astrophysical origin of this separation can¬ 
not be excluded, the authors interpret this in terms of 
DM self-interactions and infer: o-dm/uidm ~ (1.7 ± 0.7) x 
10“"^ cm^ g~^. If correct this would rule out many attractive 
DM candidates, in particular axions, neutrinos and ‘weakly 
interacting massive particles’ (WIMPs) such as supersym¬ 
metric neutralinos. 

The cross section estimated by [Massey et al.' (20151 is 
well below the weak upper bounds set by other astrophysical 
objects such as the ‘Bullet Cluster’ (IE 0657-56) which are 


typically: ajoM/muM < 1 cm^ g ^ 

(Markevitch et al. 

2004 

Randall et al.||20081 [Peter et al.j 

>012 

Rocha et al. 

2013 

Harvey et al.[[2015 

1 . [Massey et al. 

(2015 

1 argue that A3827 


is uniquely sensitive because the infall time for the galaxies 
is very long, hence the effects of DM self-interactions add up 
over a long period. They adopt a simple model for the effect 


of DM self-interactions following Williams & Saha (20111 


* felix.kahlhoefer@desy.de 


which predicts that the separation between the stars and the 
DM subhalo should grow proportional to tfnfaii thus vastly 
amplifying the effect over time. 

In this letter we argue that the model used to interpret 
the observations in terms of DM self-interactions is based 
on two questionable assumptions: 

(i) The stars and the associated DM subhalo are assumed 
to develop completely independently such that even a tiny 
difference in the acceleration they experience can lead to 
sizeable differences in their trajectories. This neglects the 
crucial fact that initially the stars are gravitationally bound 
to the DM subhalo and can only be separated from it if ex¬ 
ternal forces are comparable to the gravitational attraction 
within the system. The effect of DM self-interactions must 
therefore be at least comparable to the relevant gravitational 
effects in order to lead to an observable separation. 

(ii) The effective drag force on the DM subhalo is assumed 
to be constant throughout the evolution of the system. This 
is in contradiction with the fact that the rate of DM self¬ 
interactions typically depends both on the velocity of the 
subhalo relative to the cluster and the DM density of the 
cluster (at the position of the subhalo), both of which will 
vary along the trajectory of the subhaloj^In particular, the 


^ The assumption of a constant drag force would be justified if 
the subhalo were on a circular orbit around the centre of A3827. 
However, the observed separation between stars and DM points 
approximately in the radial direction, so this is disfavoured. 
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rate of DM self-interactions will be negligibly small as long 
as the subhalo is far away from the core of the cluster. 


These issues have previously been discussed 


by Kahlhoefer et al. (20141 (henceforth Kal4) in the 
context of merging clusters such as Abell 520, the ‘Bullet 
Cluster’ or the ‘Musket Ball Cluster’, but they apply 
equally in the present context. We apply the arguments 
from Kal4 to provide a corrected estimate of the DM 
self-interaction cross section necessary to explain the 
observations of A3827. We find values that are intriguingly 
close to existing bounds from other systems, implying that 
such systems can potentially be used to confirm or rule out 
this interpretation. 

First we provide a simple estimate of the magnitude of 
the drag force and the relevant gravitational force to extract 
the relevant self-interaction cross section. We refine this es¬ 
timate by considering a realistic trajectory for the subhalo 
and evaluating its velocity and the background DM density 
along this trajectory. We then run numerical simulations of a 
subhalo falling towards the core of A3827, accounting for the 
motion of a large number of DM test particles undergoing 
self-interactions in a time-dependent gravitational potential. 
Finally, we discuss alternatives to a simple drag force that 
are motivated from a particle physics perspective. 


2 ESTIMATES OF DM SELF-INTERACTIONS 

Let us assume that DM particles moving with velocity v 
through an ambient DM density p experience a drag force 
of the form 

Tdrag 1 d 2 /n 

- = 7 -w P . (1) 

rriDM 4 rriDM 

where a is the effective DM self-interaction cross section, 
moM is the DM mass and we have chosen the normalisation 
of a in such a way that the momentum transfer cross section 
is given by ut — 512 in analogy to the case of isotropic scat¬ 
tering (see Kal4). Such a drag force can be obtained by av¬ 
eraging over a large number of DM scattering processes with 
small scattering angle, assuming that the differential cross 
section has no velocity dependence but a strong angular de¬ 
pendence such that the probability of scattering is peaked 
in the forward direction]^ The presence of such a drag force 
will slow down any DM subhalo falling into a larger DM halo 
compared to objects experiencing only gravitational forces, 
such as the stars bound in the subhalo. 

The stars inside the subhalo, however, do not just expe¬ 
rience the gravitational attraction of the cluster but also the 
attraction of the subhalo itself. The latter contribution aims 
to reduce any separation between the subhalo and the stars 
and will therefore oppose the effect of the drag force. At 
some point, the gravitational attraction between DM sub¬ 
halo and stars will become sufficiently large to balance the 
drag force on the DM subhalo. Denoting the (average) sepa¬ 
ration between the subhalo and the stars by A, the restoring 


^ We note that such a cross section is very difficult to motivate 
from the particle physics picture. This issue will be discussed in 
more detail below. 


force will be given by 

E.h ^ GNM,h(A) 
tUstar A^ 


with Gn Newton’s constant and Msjh(A) the subhalo mass 
within radius A. As a rough estimate, we assume a constant 
density core with radius Ush = 2.7 kpc and mass Afsh = 


7 X 10^°M(T) (consistent with the mass estimate by Massey 
(2015l), so that Mah(A) = MshA^/a^j,. In order for 


et al. 


the drag force to result in a sizeable separation, we must 
require Egh/mstar < Edrag/n^DM- This inequality leads to 


5-^4 Gn Afsh A 
rriDM p 

Based on our mass model for A3827 (Appendix]^, we esti¬ 
mate p ~ 4 GeV cm“® and v ~ 1500 km s“^ at r = 15 kpc. 
Substituting these values and requiring A = 1.6 kpc yields 


-R- > 2 cm^ g \ (4) 

rriDM 

which is in some tension with the upper bound from the 
Bullet Cluster: tr/rnDM < 1.2 cm^ g”'^ (Kal4). 

If the self-interaction cross section is much smaller than 
the estimate obtained above, the stars will remain closely 
bound to the subhalo due to the overwhelming gravitational 
restoring force. In particular, it should be clear that no vis¬ 
ible separation can result from a cross section as small as 
d/moM ~ lO”"* cm^ g“^. 


2.1 One-dimensional simulations 


For the estimate above, we have assumed both the velocity 
of the subhalo v and the background density of the cluster p 
to be constant in time. In this case, the separation between 
the DM subhalo and the stars is expected to remain con¬ 
stant in time, with the equilibrium value determined by the 
condition Edrag/niDM = Esh/mstar. However as long as the 
subhalo is far away from the central region of the cluster, 
both V and p are expected to be small, so that no significant 
separation will occur. As the subhalo accelerates towards 
the central region, the separation will grow, provided that 
the drag force is always sufficiently larger than the restoring 
gravitational force. To be more realistic, we should there¬ 
fore calculate the position and velocity of the subcluster as 
a function of time. To this end, we need to also include the 
gravitational force of the cluster, acting on both the subhalo 
and the stars: 


Ecluster _ Gn Afcluster (^) 


(5) 


where Mciuater(f') is the cluster mass within radius r. As¬ 
suming a radial orbit, we can write 

.. Ecluster , Edrag 

T-sh =-h 


^sta 


mDM rriDM 
Gn -^^cluster('^sh) 


-^cluster ^ -^sh 


+ 7 r^h Pir^h) (6) 

4 TTiDM 


^star ^star 

Gn '^^cluster('^star) Gn -^^sh('^sh '^star) 




2^ 


(7) 
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Figure 1. The separation between the DM subhalo and the stars 
as a function of the subhalo distance from the centre based on a 
simple one-dimensional simulation. The vertical line indicates the 
currently observed position of the subhalo. 


where we assume the gravitational pull of the stars on the 
subhalo to be negligible. As before, we will assume a con¬ 
stant density core for the subhalo, such that Fah is pro¬ 
portional to Tah — fstar- For the cluster, however, we will 
use a more refined mass model as discussed in Appendix [X| 
and calculate the cluster mass within radius r according to 
Maiuatei{r) = dtT /J" r'^p(r') dr'. 

The set of differential equations introduced above is 
readily solved for given initial conditions. We assume that 
the subhalo starts falling towards the cluster from rest at 
an initial distance of 100 kpc and that its trajectory lies 
in the plane of the sky. It will then take approximately 
tinfaii ~ 10® yrs to reach the central region of the cluster. 
In the absence of a drag force, the subhalo will reach a ve¬ 
locity of around 1900 km s“^ at a distance of 15 kpc, which 
is similar to the velocity assumed above. 

Our results for A = Vah — rstar are shown in Fig. as 
a function of rsh taking d/mDM = 2 cm^g“^ and 3 cm^g“^. 
The observed separation at rah = 15 kpc is found to be A = 
0.8 kpc and 1.2 kpc respectively, showing that our simple 
estimate above was quite accurate. It is worth emphasising 
that the separation only becomes large once the subhalo 
comes close to the central region of the cluster. 


2.2 Three-dimensional simulations 

So far we have been using a simple model of the infalling 
subhalo with a constant density core which leads to a com¬ 
parably shallow potential and a relatively small binding en¬ 
ergy. For a more realistic halo profile, stars in the central 
region will be much more tightly bound and will therefore 
always remain very close to the peak of the DM mass dis¬ 
tribution. Indeed, the drag force on the DM subhalo can be 
interpreted as a tilt of the potential, so that loosely bound 
stars are typically more affected and will preferentially travel 
ahead of the DM subhalo (or even escape from the poten¬ 
tial) and also the minimum of the potential will be slightly 
shifted. Both effects lead to an apparent separation between 
the distribution of stars and the DM subhalo. 

Clearly, the motion of the stars in the combined poten¬ 
tial of the subhalo and the cluster cannot be modelled in 
a one-dimensional simulation. In Kal4 a three-dimensional 
simulation was developed to address this problem in the con¬ 


text of major mergers like the Bullet Cluster. Our approach 
was to treat the gravitational potential of the cluster as time- 
independent, while for the subhalo the central density and 
scale radius are allowed to vary with time and determined 
self-consistently from the simulation. Assuming an initial 
density profile, the simulation chooses a representative set 
of DM particles and stars bound to the subhalo and then 
calculates the motion of all these particles in the combined 
gravitational potential of the cluster and the subhalo. It is 
then straight-forward to add in an additional drag force af¬ 
fecting only the DM particles, based on the velocity of these 
particles and the background DM density. 

We use a Hernquist ( |1990 1 profile to model both the 
cluster and the subhalo (see Appendix]^, the advantage 
being that it has a finite central potential and the veloc¬ 
ity distribution function can be described analytically. Since 
our results are quite independent of the DM density at 
large radii, we expect very similar results for the alterna¬ 
tive Navarro-Frenk-White profile. We use the same initial 
conditions for the position and velocity of the subhalo as in 
the one-dimensional simulation. 

A significant complication of the three-dimensional sim¬ 
ulation is the need to determine the centroid of the DM sub¬ 
halo and the distribution of stars. As discussed in Kal4, it 
is inconsistent to just calculate these centroids including all 
particles which were initially bound to the subhalo, because 
doing so would include particles that have escaped from the 
DM subhalo and are now far away from the peak of the 
distribution. A realistic estimate of the separation can be 
obtained by including only particles within the projected 
iso-density contour containing 20% of the total mass of the 
original DM subhalo and ignoring regions where the sur¬ 
face density of DM particles originating from the subhalo is 
smaller than the background surface density of the cluster. 
This procedure corresponds roughly to including particles 
within 4kpc of the peak of the distribution. The same proce¬ 
dure is applied to determine the centroid of the distribution 
of stars. 

To illustrate how strongly our results depend on the 
definition of the centroid we also show the predicted sep¬ 
aration for different centroid definitions. In particular, we 
consider the case where only the innermost 5% of the total 
mass are included in the calculation of the centroid (such 
that the tails of the distribution are mostly ignored and the 
position of the centroid is determined largely by the posi¬ 
tion of the peak). While this definition may be realistic for 
the case of the stars, where the peak of the distribution can 
be accurately determined, it may not be appropriate for the 
determination of the centroid of the DM subhalo. We there¬ 
fore also consider an alternative definition, where the inner 
20% of the DM subhalo but only the inner 5% of the stars 
are included in the calculation of the centroids. 

Including the tails of the DM distribution has a fur¬ 
ther important advantage: We will see below that these tails 
can contain relevant information on the nature of DM self¬ 
interactions and should therefore not be neglected. 

Our results for the separation A are shown in Fig. as 
a function of rsh for a/m-DM = 2 cm^ g“^ and 3 cm^ g“^. 
The observed separation at rsh = 15 kpc is found to be A ~ 
1.0 kpc and ~ 1.5 kpc, in good agreement with our previous 
estimate. We also note that the predicted separation varies 
by only about 5% as we consider different definitions of the 
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Figure 2. Left: Predicted separation between the centroid of the DM subhalo and the centroid of the stars resulting from an effective 
drag force on the DM subhalo as a function of the subhalo position The solid, dashed and dot-dashed lines correspond to centroid 
calculations using the inner (20%, 20%), (5%, 5%) and (20%, 5%) of the DM subhalo and the distribution of stars, respectively (see text 
for details). Right: Histogram of the distribution of DM particles (blue) and stars (orange) as a function of radial distance r from the 
centre of the cluster at the moment in time where the radial position of the peak of the subhalo is rgh = 15 kpc (indicated by the vertical 
dashed line in the left panel). 


centroids. This indicates that not only the tails but also the 
position of the peaks differ for the two distributions. 

To confirm this expectation, we show in the right panel 
of Fig. the distribution of the stars and the DM parti¬ 
cles along the radial direction. Indeed, one can see that the 
peaks of the two distributions are slightly shifted. Further¬ 
more the tail of the distribution of stars is enhanced in the 
forward direction due to stars that have escaped from the 
gravitational potential of the subhalo. Most of the remain¬ 
ing stars, however, remain bound to the subhalo and would 
return to their equilibrium position if the subhalo survived 
the passage through the central region of the cluster. 


2.3 Alternative interpretations 


We briefly discuss alternative particle physics models for DM 
self-interactions and their interpretation in the present con¬ 
text. The two most natural forms of DM self-interactions are 
long-range interactions, which can arise for example from a 


massless mediator (Feng et al. 20091, and contact interac¬ 


tions. Long-range interactions typically lead to a drag force 
of the form 


- drag 


moM 4 moM 


( 8 ) 


with a/moM < 10“^^ cm^ g“^ from the observation that 
Milky Way satellites have survived up to the present 
day ( [Gnedin fc Qstriker||200T| . Using an estimate similar 
to the one discussed above, one can immediately see that 
for the DM densities and the velocities under consideration 
such a drag force cannot give an observable effect. More¬ 
over, the inverse velocity dependence implies that DM self¬ 
interactions would actually be suppressed as the subhalo 
approaches the central region of the cluster, in contrast to 
what is required in order to obtain large effects from DM 
self-interactions. However a massless mediator can induce 


collisionless shocks, which Heikinheimo et al. (20151 argue 


can potentially explain the features observed in A3827. 

For contact interactions, an effective description in 


terms of a drag force is not possible, because each indi¬ 
vidual DM particle will only experience a small number of 
collisions (if any). The reason is that in each collision the 
momentum transfer is so large that these collisions must be 
rare to avoid observational bounds from the survival of sub¬ 
haloes and halo ellipticity ([Peter et al.||201^ [Rocha et al.[ 


20131. Nevertheless, a separation between the DM subhalo 


and stars can in principle also occur in the case of contact 
interactions. In this case, the separation is not due to stars 
leaving the DM subhalo in the forward direction, but due to 
DM particles receiving large momentum transfer and leaving 
the DM subhalo in the backward direction. In other words, 
rare DM self-interactions can affect a tail of DM particles, 
which shifts the centroid of the DM distribution relative to 
the distribution of stars, while the peak of the DM distri¬ 
bution, which is dominated by DM particles that have not 
experienced any self-interactions, remains coincident with 
the peak of the distribution of stars. 


This scenario can also be investigated using the simula¬ 
tions discussed above (for details, see Kal4). Our results are 
shown in Fig. [^for cr/rriDM = l.Ocm^ g“^ and 1.5cm^ g“^. 
For larger self-interaction cross sections, the evaporation 
rate of the DM subhalo becomes so large that it is unlikely 
to have survived up to its present position. Including the 
inner 20% of the DM subhalo and the stars in the centroid 
calculation, the observed separations at rah = 15 kpc are 
found to be A = 1.0 kpc and 1.6 kpc, respectively. As in 
the case of a drag force, there is tension between the ob¬ 
servation of a separation in A3827 and other constraints on 
the DM self-interaction cross section, which are typically 


o-DM/moM < 1 cm^ g ^ ( 

Markevitch et al. 

2004 

Randall 

et al.||2008| [Peter et al.||2012| [Rocha et al. 

[20131 [Harvey 


et al.|2015[ ). 


In contrast to the case of an effective drag force we find 
that the predicted separation depends sensitively on the def¬ 
inition of the centroids. Including only the innermost 5% of 
the subhalo and the stars (in order to determine the approx¬ 
imate position of the respective peaks) yields a somewhat 
smaller predicted separation of 0.8 kpc and 1.3 kpc. This 
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Figure 3. Same as Fig.j^but for the case of rare self-interactions. In this case, the DM subhalo develops a tail in the backward direction, 
while the distribution of stars remains largely unchanged. The peaks of the two distributions remain largely coincident and therefore the 
separation depends sensitively on whether the tails of the distributions are included in the centroid calculation. 


observation indicates that the separation is mainly due to 
differences in the shapes of the two respective distributions, 
while the peaks of the distributions remain coincident. This 
expectation is confirmed by the histrograms shown in the 
right panel of Fig. 

An important conclusion is that the case of contact in¬ 
teractions can potentially be distinguished from the case of 
an effective drag force by studying in detail the shape of 
the DM subhalo and the relative position of the peaks of 
the two distributions. In the case of contact interactions, the 
DM subhalo is expected to be deformed due to the scattered 
DM particles leaving the subhalo in the backward direction, 
such that the position of the centroid depends sensitively 
on which particles are included in the calculation. For an 
effective drag force, on the other hand, we expect the DM 
subhalo to retain its shape, while the distribution of stars 
will be both shifted and deformed. 


3 DISCUSSION 

In this letter we have discussed a possible interpretation in 
terms of DM self-interactions of an observed separation of 
about 1.6kpc between the DM halo of a galaxy and its stars. 
Using several increasingly refined methods we estimate that 
the self-interaction cross sections necessary to explain this 
effect are of order a/m-DM ~ 3 cm^ g“^ for the case of an ef¬ 
fective drag force (proportional to the square of the velocity) 
or a/moM ~ 1.5cm^g“^ for the case of contact interactions. 
Both of these values are in tension with the upper bounds on 
DM self-interactions from other astrophysical observations. 

We have modelled the system under consideration 
rather simply, but our derived estimates are conservative 
since a more refined analysis (e.g. considering other possi¬ 
ble trajectories for the subhalo) would predict a higher self¬ 
interaction cross section for the same observed separation. It 
should therefore be clear that A3827 is no more sensitive to 
DM self-interactions than other systems considered in this 
context and can certainly not be used to probe cross-sections 
as small as aoM/moM ~ 10“^ cm^ g“^. Further studies of 
such systems are imperative to establish if the indication 
from A3827 for a non-zero self-interaction cross section (of 
order 1 cm^ g~^) is indeed correct. 
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APPENDIX A: A CRUDE MODEL OF A3827 


For simplicity, we model both the central region of clus¬ 
ter A3827 and the galactic subhalo called N1 b y |Williams 
& Saha (20111 and Massey et al. (20151 using a Hernquist 


(19901 profile: 


p(r) = 


M 


2n r {a + r)® 


(Al) 
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For the cluster, we determine M^uster and Uduster by fitting 
the mass dis tribution to the one given in figure B2 of |Masse3^ 
(20151. This procedure leads to Mduster = 7x 


et al. 


and Ociuster = 60 kpc. As suggested in Massey et al. (20151 


we assume a constant-density core with radius 8 kpc. We 
have checked that this model reproduces the observed pro¬ 
jected mass in the central region and the observed velocity 
dispersion with reasonable accuracy. 

For the subhalo, we find that taking M^h = 5 x 10^^ Mq 
and Osh = 7 kpc enables us to reproduce satisfactorily the 
observed velocity dispersion and the projected mass in the 
central region of the subhalo. 
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